Lorentz group in ray optics 



Sibel Ba§kal ^ 

Department of Physics, Middle East Technical University, 
06531 Ankara, Turkey 

Elena Georgieva^ 

Science Systems and Applications, Inc., Lanham, MD 20771, U.S.A., and 
National Aeronautics and Space Administration, 
Goddard Space Flight Center, Laser and Electro-Optics Branch, 
Code 554, Greenbelt, Maryland 20771, U.S.A. 

Y. S. Kim^ 

Department of Physics, University of Maryland, 
College Park, Maryland 20742, U.S.A. 

Marilyn E. Noz ^ 

Department of Radiology, New York University, New York, New York 

10016, U.S.A. 

Abstract 

It has been almost one hundred years since Einstein formulated 
his special theory of relativity in 1905. He showed that the basic 
space-time symmetry is dictated by the Lorentz group. It is shown 
that this group of Lorentz transformations is not only applicable to 
special relativity, but also constitutes the scientihc language for op- 
tical sciences. It is noted that coherent and squeezed states of light 
are representations of the Lorentz group. The Lorentz group is also 
the basic underlying language for classical ray optics, including polar- 
ization optics, interferometers, the Poincaree sphere, one-lens optics, 
multi-lens optics, laser cavities, as well multilayer optics. 



-'^electronic address:baskal@newton. physics. metu.edu.tr 
^electronic address: egeorgie@pop500.gsfc.nasa.gov 
•^electronic address: yskini@physics.umd.edu 
^electronic address: noz@nucmed.med.nyu.edu 



1 



1 Introduction 

Since Einstein formulated special relativity in 1905, the basic space-time 
symmetry has been that of the Lorentz group. He established the energy- 
momentum relation which is valid for slow massive particles, high-speed mas- 
sive particles, and massless particles. Einstein formulated this relation ini- 
tially for particles without internal space-time structures, but it is widely 
accepted that it is valid for all particles including particles with spin and/or 
internal space-time extension. 

In 1939, Wigner published his most fundamental paper dealing with inter- 
nal space-time symmetries of relativistic particles PP. In this paper, Wigner 
introduced the Lorentz group to physics. Furthermore, by introducing his 
"little groups," Wigner provided the framework for studying the internal 
space-time symmetries of relativistic particles. The scientific contents of this 
paper have not yet been fully recognized by the physics community. We are 
writing this report as a continuation of the work Wigner initiated in this 
history-making paper. 

While particle physicists are still struggling to understand internal space- 
time symmetries of elementary particles, Wigner's Lorentz group is becoming 
useful to many other branches of physics. Among them are optical sciences, 
both quantum and classical. In quantum optics, the coherent and squeezed 
states are representations of the Lorentz group Recently, the Lorentz 
group is becoming the fundamental language for classical ray optics. It is 
gratifying to note that optical components, such as lenses, polarizers, in- 
terferometers, lasers, and multi-layers can all be formulated in terms of the 
Lorentz group which Wigner formulated in his 1939 paper. Classical ray 
optics is of course a very old subject, but we cannot do new physics without 
measurements using optical instruments. Indeed, classical ray optics consti- 
tutes Wigner's frontier in physics. 

The word "group theory" sounds like abstract mathematics, but it is 
gratifying to note that Wigner's little groups can be formulated in terms of 
two-by-two matrices, while classical ray optics is largely a physics of two-by- 
two matrices. The mathematical correspondence is straight-forward. 

In order to see this point clearly, let us start with the following classic 
example. The second-order differential equation 
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Table 1: Further contents of Einstein's E = mc^. Massive and massless par- 
ticles have different energy-momentum relations. Einstein's special relativity 
gives one relation for both. Wigner's httle group unifies the internal space- 
time symmetries for massive and massless particles. The quark model and 
the parton model can also be combined into one covariant model. 
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is applicable to a driven harmonic oscillator with dissipation. This can also be 
used for studying an electronic circuits consisting of inductance, resistance, 
capacitance, and an alternator. Thus, it is possible to study the oscillator 
system using an electronic circuit. Likewise, an algebra of two-by-two ma- 
trices can serve as the scientific language for several different branches of 
physics, including special relativity, ray optics, and quantum optics. 

There are many physical systems which can be formulated in terms of 
two-by-two matrices. If we restrict that their determinant be one, there is a 
well established mathematical discipline called the group theory of SL{2, C). 
This aspect was noted in the study of Lorentz transformations. In group 
theoretical terminology, the group SL{2, C) is the universal covering group 
of the group of the Lorentz group. In practical terms, to each two-by-two 
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matrix, there corresponds one four-by-four matrix which performs a Lorentz 
transformation in the four-dimensional Minkowskian space. Thus, if a phys- 
ical system can be explained in terms of two-by-two matrices, it can be 
explained with the language of Lorentz transformations. Furthermore, the 
physical system based on two-by-two matrices can serve as an analogue com- 
puter for Lorentz transformations. 

Optical filters, polarizers, and interferometers deal with two independent 
optical rays. They superpose the beams, change the relative phase shift, 
and change relative amplitudes. The basic language here is called the Jones 
matrix formalism, consisting of the two-by-two matrix representation of the 
SL{2,C) group OH]. The four-by-four Mueller matrices are derivable from 
the two- by-two matrices of SL{2,C). 

For these two-beam systems, the Poincare sphere serves as an effective 
language jSJIHl. Since the two-beam system is described by the Lorentz group, 
the Poincare sphere is necessarily a representation of the Lorentz group. We 
shall use this sphere to study the degree of coherence between the beams. 

Para-axial lens optics can also be formulated in terms of two-by-two ma- 
trices, applicable to the two-component vector space consisting of the dis- 
tance from the optical axis and the slope with respect to the axis. The lens 
and translation matrices are triangular, but they are basically representations 
of the Sp{2) group which is the real subgroup of the group SL{2, C) [3 |H]. 

Laser optics is basically multi-lens lens optics. However, the problem 
here is how to get a simple mathematical expression for the system of a large 
number of the same lens separateled by equal distance. Here again, group 
theory simplifies calculations p. 

In multi-layer optics, we deal with two optical rays moving in opposite 
directions. The standard language in this case is the S-matrix formalism [TUj . 
This is also a two-by-two matrix formalism. As in the case of laser cavities, 
the problem is the multiplication of a large number of matrix chains |TT1 IT^ . 

It is shown in this report that the two-by-two representation of the six- 
parameter Lorentz group is the underlying common scientific language for all 
of the instruments mentioned above. While the abstract group theoretical 
ideas make two-by-two matrix calculations more systematic and transparent 
in optics, optical instruments can act as analogue computers for Lorentz 
transformations in special relativity. It is gratifying to note that special 
relativity and ray optics can be formulated as the physics of two-by-two 
matrices. 

In Sec. 121 we explain how the Lorentz group can be formulated in terms 
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of four-by-four matrices. It is shown that the group can have six independent 
parameters. In Sec. |21 we explain how it is possible to formulate the Lorentz 
group in terms of two-by-two matrices. It is shown that the four-by-four 
transformation matrices can be constructed from those those of the two-by- 
two matrices. 

In Sec. m we discuss the historical significance of Wigner's 1939 paper [T] 
on the Lorentz group and its application to the internal space-time symme- 
tries of relativistic particles. In Sec. we present the basic building blocks for 
the two-by-two representation of the Lorentz group in terms of the matrices 
commonly seen in ray optics. 

In Sees. and IHl we study the polarizations for the two-beam system. It 
is shown that both the Jones and Mueller matrices are representations of the 
Lorentz group. The role of the Stokes parameters is discussed in detail. It is 
shown in Sec. IHlthat the Poincare sphere is a graphical representation of the 
Poincare group, and serves as a device to describe the degree of coherence. 

In Sees, [ini im [m andlT^ we discuss one-lens system, multi-lens system, 
laser cavities, and multi-layer optics, respectively. In all of these sections, 
the central scientific language is the Lorentz group. 

In Appendix El we expand the content of the third row of Table ^ It is 
noted that the covariant harmonic oscillator formalism can unify the quark 
model for slow hadrons with the parton model for ultra-relativistic quarks. It 
is then shown that the same oscillator formalism serves as the basic scientific 
language for squeezed states of light. 

In Appendix ^ it is shown that the Lie-group method, in terms of the 
generators, is not the only method in constructing group representations. 
For the rotation group and the three-parameter subgroups of the Lorentz 
group, it is simpler to start with the minimum number of starter matrices. 
For instance, while there are three generators for the rotation group in the 
Lie approach, we can construct the most general form of the rotation matrix 
from rotations around two directions, as Goldstein constructed the Euler 
angles jlSl- 

In Appendix O it is noted that the four- by-four matrices are real, their 
two-by-two counterparts are complex. However, there is a three-parameter 
real subgroup called 5*^(2). It is shown that the complex subgroup SU{1, 1) 
is equivalent to 5*^(2) through conjugate transformation. 
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2 Lorentz Transformations 



Let us consider the space-time coordinates {ct, z,x,y). Then the rotation 
around the z axis is performed by the four-by-four matrix 
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Likewise, we can write down the generators of rotations Ji and J2 around 
the X and y axes respectively. 
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These three generators satisfy the closed set of commutations relations 

[Ji, Jj] = iCijkJk- (5) 

This set of commutation relations is for the three-dimensional rotation group. 
The Lorentz boost along the z axis takes the form 
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Likewise, we can write generators of boosts Ki and K2 along the x and y 
axes respectively, and they take the form 



/O i 0\ 
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Vo 0/ 



fO i\ 





\i OJ 



(8) 



These boost generators satisfy the commutation relations 

[Ji,Kj] = ieijkKk, [Ki,Kj] = -ieijkJk- 



(9) 



Indeed, the three rotation generators and the three boost generators sat- 
isfy the closed set of commutation relations given in Eq.(0) and Eq.Q. These 
three commutation relations form the starting point of the Lorentz group. 
The generators given in this section are four-by-four matrices, but they are 
not the only set satisfying the commutation relations. We can construct also 
six two-by-two matrices satisfying the same set of commutation relations. 
The group of transformations constructed from these two-by-two matrices 
is often called SL{2, c) or the two-dimensional representation of the Lorentz 
group. Throughout the present paper, we used the two-by-two transforma- 
tion matrices constructed from the generators of the SL{2, c) group. 



3 Spinors and Four-vectors 

In Sec. m we started with four-by-four transformation matrices applicable 
to the four-dimensional space-time. We then ended up with a set of closed 
commutation relations for the six generators consisting of three rotation and 
three boost generators. These generators are in the form of four-by-four 
matrices. 

In this section, we shall see first that there is a set of two-by-two ma- 
trices satisfying the same set of commutation relations, constituting the 
two-by-two representations of the Lorentz group. The representation so 
constructed is called the SL{2, c) group, or the universal covering group 
of the Lorentz group. The transformation matrices are applicable to two- 
component SL{2, C) spinors. The algebraic property of this two-by-two rep- 
resentation is the same as that of the four-by-four representation. 

As we shall see in this paper, the spinors and the four-vectors correspond 
to Jones vectors and Stokes parameters respectively in polarization optics |3] . 
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The question then is whether we can construct the four-vector from the 
spinors. In the language of polarization optics, the question is whether it is 
possible to construct the coherency matrix IHj from the Jones vector. 

With this point in mind, let us start from the following form of the Pauli 
spin matrices: 

These matrices are written in a different convention. Here is imaginary, 
while (72 is imaginary in the traditional notation. Also in this convention, we 
can construct three rotation generators 

Jt = ^CTi, (11) 
which satisfy the closed set of commutation relations 

[Ji, Jj] = itijkJk- (12) 
We can also construct three boost generators 

= (13) 
which satisfy the commutation relations 

[Ki,Kj\ = -ieijkJk- (14) 

The Ki matrices alone do not form a closed set of commutation relations, 
and the rotation generators Jj are needed to form a closed set: 

[Ji,Kj] = ieijkKk. (15) 

The six matrices Ji and Ki form a closed set of commutation relations, 
and they are like the generators of the Lorentz group applicable to the (3 
+ l)-dimensional Minkowski space. The group generated by the above six 
matrices is called SL{2, c) consisting of all two-by-two complex matrices with 
unit determinant. 

Let us write the two-by-two transformation matrix as 



a P 
7 6 



(16) 
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This matrix has four complex elements with eight real parameters. However, 
the six generators are all traceless and, the determinant of the matrix has to 
be one. Thus, it has six independent parameters. 

In order to construct four-vectors, we need two different spinor represen- 
tations of the Lorentz group. Let us go to the commutation relations for 
the generators given in Eqs. (jl2|) . (jl4j) and (jl5p . These commutators are not 
invariant under the sign change of the rotation generators Jj, but are invari- 
ant under the sign change of the squeeze operators Ki. Thus, to each spinor 
representation, there is another representation with the squeeze generators 
with opposite sign. This allows us to construct another representation with 
the generators: 

1 • —i 

Ji = -o-i, Ki = —Gi. (17) 

We call this representation the "dotted" representation. If we write the 
transformation matrix L of Eq. (jl6|) in terms of the generators as 




then the transformation matrix in the dotted representation becomes 




In both of the above matrices, Hermitian conjugation changes the direction of 
rotation. However, it does not change the direction of boosts. We can achieve 
this only by changing L to L, and we shall call this the "dot" conjugation. 
Likewise, there are two different set of spinors. Let us use write 




for the spinor in the undotted and dotted representations respectively. Then 
the four- vectors are constructed as [HI [T3] 

uu = — (x — iy\ f -y = (x + iy), 

uv = {t^z), vu = -{t-z) (21) 

leading to the matrix 

V=f"^ -"^') = fM(^ -^), (22) 
V vv —vu J \v J 
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where u and ii are one if the spin is up, and are zero if the spin is down, 
while V and v are zero and one for the spin-up and spin-down cases. 

If the two-by-two matrix of Eq. is apphcable to the column vector 
of Eq. ()20|) . what is the transformation matrix applicable to the row vector 
{v, —u) from the right-hand side? It is the transpose of the matrix applicable 
to the column vector iv, — u). We can obtain this column vector from 



V 

—u 



(23) 



by applying to it the matrix 



9 = -^^3={l q)- (24) 

This matrix also has the property 

gcTig" = - ((Ti) , (25) 

where the superscript T means the transpose of the matrix. The transforma- 
tion matrix applicable to the column vector of Eq. lj^Hj) is L of Eq. p9p . Thus 
the matrix applicable to the row vector [v, — u) in Eq.()22|) is 

{g-^Lgf = g-^L'g. (26) 

This is precisely the Hermitian conjugate of L. 
Let us now write the V matrix of Eq. (j^^ as 



,x + iy t — z 

where the set of variables {x, y, z, t) is transformed like a four- vector under 
Lorentz transformations. Then the Lorentz transformation on V can be 
performed as 

V = LVL\ (28) 

where the transformation matrix L is that of Eq.()16p. As we have seen in this 
section, the construction of four-vectors from the two-component spinors is 
not a trivial task, but has been discussed in the literature Likewise, 
it is possible to construct four-by-four Lorentz transformation matrices from 
the two- by-two matrix of Eq. (fTBj) ITT] . 

As we shall see in the present paper, the matrix of the form Eq. fl22j) or 
Eq. appears in optics as the coherency matrix and the density matrix for 
two beam systems such as polarization optics and interferometers. 
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4 Wigner's Little Groups 



In his 1939 paper, Wigner introduced his httle groups in order to study 
the internal space-time structures of relativistic particles. The little group 
is the maximal subgroup of the Lorentz group which leaves the momentum 
of the particle invariant. The little groups, originally developed for studying 
relativistic symmetries of particles, are now becoming an important scientific 
language for ray optics as we shall see in this paper. 

Let us first discuss the physical basis of Wigner's little groups. If the 
speed of a particle is much smaller than that of light, energy-momentum 
relation is E = p^/2mo + moc^. If the speed is close to that of light, the 
relation is E = cp. These two different relations can be combined into one 
covariant formula E"^ = m^c^ + p^c^. This aspect of Einstein's E = mc^ is 
well known, as indicated in Table ^ 

In addition, particles have internal space-time variables. Massive particles 
have spins while massless particles have their helicities and gauge variables. 
Our first question is whether this aspect of space-time variables can be unified 
into one covariant concept. The answer to this question is Yes. Wigner's little 
group does the job, also as indicated Table [TJ 

particles can also have space-time extensions. For instance, in the quark 
model, hadrons are bound states of quarks. However, the hadron appears 
as a collection of partons when it moves with speed close to the velocity of 
light. Quarks and partons seem to have quite distinct properties. Are they 
different manifestations of a single covariant entity? This is one of the most 
pressing issues in high-energy particles physics. The third row of Table ^ 
addresses this question. 

The mathematical framework of this program was developed by Eugene 
Wigner in 1939 ilj. He constructed the maximal subgroups of the Lorentz 
group whose transformations will leave the four-momentum of a given particle 
invariant. These groups are known as Wigner's little groups. Thus, the 
transformations of the little groups change the internal space-time variables 
of the particle, while leaving its momentum invariant. 

In his paper jTj, Wigner shows that, for each massive particle, there is a 
Lorentz frame in which the particle is at rest. Then the three-dimensional 
rotation group leaves its momentum invariant, while changing the direction 
of its spin. Thus, the httle group for this particular case is 0(3). In other 
Lorentz frames, the little group is a Lorentz-boosted 0(3) group. 

For a massless particle, Wigner notes that there are no Lorentz frames in 
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which the particle is at rest. The best we can do is to ahgn its momentum to 
a given axis, or the z axis. Then, its momentum is invariant under rotations 
around the z axis. In addition, Wigner found two more transformations which 
leave the momentum invariant. The physics of these additional degrees was 
not explained in Wigner's original paper , but their mathematics has been 
worked out. The physics of these degrees of freedom was later determined to 
be that of gauge transformations (THj, and its complete understanding was 
not achieved until 1990 [T^ . 

The group of Lorentz transformations consists of three boosts and three 
rotations. The rotations therefore constitute a subgroup of the Lorentz 
group. If a massive particle is at rest, its four-momentum is invariant under 
rotations. Thus the little group for a massive particle at rest is the three- 
dimensional rotation group. Then what is affected by the rotation? The 
answer to this question is very simple. The particle in general has its spin. 
The spin orientation is going to be affected by the rotation! If we use the 
four- vector coordinate (ct, z^ x, y), the four-momentum vector for the particle 
at rest is (moc, 0,0,0), and the three-dimensional rotation group leaves this 
four- momentum invariant. This little group is generated by the three rota- 
tion generators given in Eq.(j2I) and Eq.(0]). They satisfy the commutation 
relations for the rotation group given Eq.(j3)). 

If the rest-particle is boosted along the z direction, it up a non-zero 
momentum component along the same direction. The 0(3) generators will 
also be boosted. The boost matrix takes the form of Eq.(jni), and the boosted 
generators will be 

J[ = B Ji B-\ (29) 

and this boost will not change the commutation relation of Eq.(0). 

For a massless particle moving along the Z direction, Wigner observed 
that the little group generated by the rotation generator around the z axis, 
namely J3 of Eq.ijS}, and two other generators which take the form 
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(30) 



If we use Ki for the boost generator along the i-th axis, these matrices can 
be written as 

N, = K,-J2, N2 = K2 + Ji, (31) 
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with Ki and K2 given in Eq.®. 

The generators J3, A^^i and N2 satisfy the following set of commutation 
relations. 

[iVi, N2] = 0, [J3, iVi] = ^N2, [J3, N2] = -tm. (32) 

In order to understand the mathematical basis of the above commutation 
relations, let us consider transformations on a two-dimensional plane with 
the xy coordinate system. We can then make rotations around the origin 
and translations along the x and y directions. If we write these generators 
as L, and Py respectively, they satisfy the commutation relations ^7j 

[P., Py] = 0, [L, P,] = iPy, [L, Py] = -iP^. (33) 

This is a closed set of commutation relations for the generators of the E{2) 
group, or the two-dimensional Euclidean group. If we replace Ni and A''2 
of Eq. ()32|) by Px and Py, and J3 by L, the commutations relations for the 
generators of the E{2)-\\ke little group becomes those for the £'(2)-like little 
group. This is precisely why we say that the little group for massless particles 
are like E{2). 

It is not difficult to associate the rotation generator J3 with the helicity 
degree of freedom of the massless particle. Then what physical variable is 
associated with the Ni and N2 generators? Indeed, Wigner was the one who 
discovered the existence of these generators, but did not give any physical 
interpretation to these translation-like generators in his original paper p. 
For this reason, for many years, only those representations with the zero- 
eigenvalues of the operators were thought to be physically meaningful 
representations [20] • It was not until 1971 when Janner and Janssen reported 
that the transformations generated by these operators are gauge transforma- 
tions ^H|- The role of this translation-like transformation has also been 
studied for spin-1/2 particles, and it was concluded that the polarization of 
neutrinos is due to gauge invariance PT| . 

The 0(3)-like little group remains 0(3)-like when the particle is Lorentz- 
boosted. Then, what happens when the particle speed becomes the speed of 
light? The energy-momentum relation E"^ = rn^c^ + p^c^ become E = pc. 
Is there then a limiting case of the 0(3)-like little group? Since those little 
groups are like the three-dimensional rotation group and the two-dimensional 
Euclidean group respectively, we are first interested in whether E{2) can be 
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obtained from 0(3). This will then give a clue as to how to obtain the E{2)- 
like little group as a limiting case of 0(3)-like little group. With this point 
in mind, let us look into this geometrical problem. 

In 1953, Inonii and Wigner formulated this problem as the contraction 
of 0(3) to E{2) [22] • Let us see what they did. We always associate the 
three-dimensional rotation group with a spherical surface. Let us consider a 
circular area of radius one kilometer centered on the north pole of the earth. 
Since the radius of the earth is more than 6,450 times longer, the circular 
region appears to be fiat. Thus, within this region, we use the E{2) symmetry 
group. The validity of this approximation depends on the ratio of the two 
radii. 

How about then the little groups which are isomorphic to 0(3) and E{2)1 
It is reasonable to expect that the E{2)-\ike little group can be obtained as 
a limiting case of the 0(3)-like little group for massless particles. In 1981, it 
was observed by Bacry and Chang [23] and by Ferrara and Savoy [21| that 
this limiting process is the Lorentz boost to infinite- momentum frame. 

In 1983, it was noted by Han et al that the large-radius limit in the 
contraction of 0(3) to E{2) corresponds to the infinite-momentum limit for 
the case of the 0(3)-like little group to the i?(2)-like little group for mass- 
less particles. They showed that transverse rotation generators become the 
generators of gauge transformations in the limit of infinite momentum |25j . 

Let us see how this happens. The J3 operator of Eq.Q, which generates 
rotations around the z axis, is not affected by the boost conjugation of the 
B matrix of Eq. (f^^ . On the other hand, the Ji and J2 matrices become 

A^i = lim e-''B-^J2B, N2 = lim -e^'^B'^JiB, (34) 

r]—*oo 77— >oo 

which are given in Eq. (j3U|) . The generators A^^i and N2 are the contracted 
J2 and Ji respectively in the infinite-momentum. In 1987, Kim and Wigner 
studied this problem in more detail and showed that the little group for 
massless particles is the cylindrical group which is isomorphic to the -£"(2) 
group [2n]- 

This completes the second row in Table ^ where Wigner's little group 
unifies the internal space-time symmetries of massive and massless particles. 
The transverse components of the rotation generators become generators of 
gauge transformations in the infinite- momentum limit. 

Let us go back to Table I given in Seed As for the third row for relativis- 
tic extended particles, the most efficient approach is to construct represen- 
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tations of the little groups using the wave functions which can be Lorentz- 
boosted. This means that we have to construct wave functions which are 
consistent with all known rules of quantum mechanics and special relativity. 
It is possible to construct harmonic oscillator wave functions which satisfy 
these conditions. We can then take the low-speed and high-speed hmits of 
the covariant harmonic oscillator wave functions for the quark model and 
the parton model respectively. This aspect was extensively discussed in the 
literature IIZ1I2ZI, and is beyond the scope of the present report. 

However, it is important to note that the covariant harmonic oscillator 
formalism use for this purpose can serve as the Fock space description for 
the squeezed state of light |2H], we give a brief discussion of this aspect in 
Appendix |X[ entitled "Squeezed Harmonic Oscillators." 

In this section, we discussed Wigner's little groups applicable to internal 
space-time symmetries of relativistic particles. However, as we shall see in 
this paper, these little groups play an important role in understanding ray 
optics. Conversely, optical configurations in ray optics can serve as analogue 
computers for space-time symmetries in particle physics. 



5 Polarization Optics 



Let us consider two optical beams propagating along the z axis. We are then 
led to the column vector: 



_ f Ai exp {—i {kz — out + (f)i)) 
'ip2{z,t)J \ A2 exp {-i (kz - Lut + (1)2)) 



(35) 



We can then achieve a phase shift between the beams by applying the two 
by-two matrix: 

e"*'^/^ 



(36) 

If we are interested in mixing up the two beams, we can apply 

(37) 



cos(e/2) -sin(^/2)\ 

sin(^/2) cos(e/2) ) 

to the column vector. 

If the amplitudes become changed either by attenuation or by reflection, 
we can use the matrix 

•e'?/2 

e~^/^ 



(38) 
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for the change. In this paper, we are deahng only with the relative ampli- 
tudes, or the ratio of the amplitudes. As we shall see in Sec. IHl the above 
two-by-two matrices have their corresponding four-by-four matrices respec- 
tively. 

Repeated applications of these matrices lead to the form 



where the elements are in general complex numbers. The determinant of this 
matrix is one. Thus, the matrix can have six independent parameters. 

This matrix takes the identical form as the two-by-two matrix given in 
Eq.(jnni)- However, the construction process is different. The Lie-group gen- 
erators are used for Eq. ljH^ . while Eq.(|H^ is constructed from repeated ap- 
plications of the transformation matrices. The difference between these two 
methods is discussed in Appendix IHl 

In either case, this matrix is the most general form of the matrices in the 
SL{2, c) group, which is known to be the universal covering group for the 
six-parameter Lorentz group. This means that, to each two-by-two matrix of 
SL{2,c), there corresponds one four-by-four matrix of the group of Lorentz 
transformations applicable to the four-dimensional Minkowski space [17 . It 
is possible to construct explicitly the four-by-four Lorentz transformation 
matrix from the parameters a, /3, 7, and 6. This expression is available in the 
literature [TTj, and we consider here only special cases. 

The four-by-four representation of the Lorentz group can be constructed 
from the two-by-two representation of the SL{2, c) group, which is known as 
the universal covering group of the six-parameter Lorentz group. This aspect 
is discussed in Sec. El 

Let us go back to Eq. ljH^ . the SL{2, c) group represented by this matrix 
has many interesting subgroups. If the matrices are to be Hermitian, then the 
subgroup is SU{2) corresponding to the three-dimensional rotation group. If 
all the elements are real numbers, the group becomes the three-parameter 
Sp{2) group. This subgroup is equivalent to SU{1, 1) which is the primary 
scientific language for squeezed states of light [3 12H] • 

We can also consider the matrix of Eq. ()39p when one of its off-diagonal 
elements vanishes. Then, it takes the form 




(39) 




(40) 
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where 7 is a complex number with two real parameters. In 1939 PP, Wigner 
observed this form as one of the subgroups of the Lorentz group. He observed 
further that this group is isomorphic to the two-dimensional Euclidean group, 
and that its four-by-four equivalent can explain the internal space-time sym- 
metries of massless particles including the photons. 

In ray optics, we often have to deal with this type of triangular matrices, 
particularly in lens optics and stability problems in laser and multi-layer 
optics. In the language of mathematics, dealing with this form is called the 
Iwasawa decomposition [U 121] • 

If the Jones matrix contains all the parameters for the polarized light 
beam, why do we need the mathematics in the four- dimensional space? The 
answer to this question is well known. In addition to the basic parameter 
given by the Jones vector, the Stokes parameters give the degree of coherence 
between the two rays. 

Let us go back to the Jones spinor of Eq. , and construct the quantities: 

Su =< Ipllpl >, 5*12 =< 1plij2 >, 

S21 =< 1p*2lpl >, ^22 =< ^2*^2 > • (41) 

Then the Stokes vector consists of 

'S'o = Sii + 5*22, >S'i = 5*11 — 5*22, 

^2 = ^12 + ^21, S3 = -^{Su + S22)■ (42) 

The four-component vector 

(^0,^1,^2,^3) (43) 

transforms like the space-time four-vector (ct, z, x, y) under Lorentz trans- 
formations. The Mueller matrix is therefore like the Lorentz-transformation 
matrix. 

As in the case of special relativity, let us consider the quantity 

= - Si - Si - Si (44) 

Then M is like the mass of the particle while the Stokes four-vector is like 
the four-momentum. 
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If M = 0, the two-beams are in a pure state. As M increases, the system 
becomes mixed, and the entropy increases. If it reaches the value of 5*0, 
the system becomes completely random. It is gratifying to note that this 
mechanism can be formulated in terms of the four-momentum in particle 
physics jl]. 

Although we can borrow all the elegant mathematical identities of the 
two-by-two representations of the Lorentz group, this formalism does not 
allow us to describe the loss of coherence within the interferometer system. 
In order to study this effect, we have to construct the coherency matrix: 

C={^'' o'^)- (45) 

\021 022 / 

Under the optical transformations discussed in this section, the coherency 
matrix is transformed as 

C = GC G\ (46) 

as in the case of the Lorentz transformation given in Eq. ()28|) . 

Using this formalisms based on the Lorentz group, we can discuss the 
group theoretical property of polarization optics in detail ^j. However, po- 
larization optics is based on two independent beams propagating in the same 
direction. Since interferometers are also based on the same optical system, 
we can continue our discussion in the following section on interferometers. 



6 Interferometers 

Typically, one beam is divided into two by a beam splitter. We can write 
the incoming beam as 

Then, the beam splitter can be written in the form of a rotation matrix jST] : 

/cos(^^/2) -sin(^^/2)\ 
^^^^ - I sin(^^/2) cos(^/2) ) ' ^^^^ 



which transforms the column vector of Eq. ()47|) into 

,7-U \ 

(49) 



- ( [cos(^/2)]exp{i(A;2 -cut)} 
i^2J ~ V-[sin(e/2)]exp{i(A;z -cut)} 
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The first beam ipi of Eq. (|47j) is now split into ipi and ip2 of Eg . (1491) . The 
intensity is conserved. If the rotation angle 6 is -7r/4, the initial beam is 
divided into two beams of the same intensity and the same phase [32j. 

These two beams go through two different optical path lengths, resulting 
in a phase difference. If the phase difference is 0, the phase shift matrix is 

W)=('o (50) 

When reflected from mirrors, or while going through beam splitters, there 
are intensity losses for both beams. The rate of loss is not the same for the 
beams. This results in the attenuation matrix of the form 



e-"^ J-"" [ e-^'/2 ' ^^^^ 

with rj = rj2 — r]i . This attenuator matrix tells us that the electric fields 
are attenuated at two different rates. The exponential factor e"*^''^"'"''^)/^ re- 
duces both components at the same rate and does not affect the degree of 
polarization. The effect of polarization is solely determined by the squeeze 
matrix 

In the detector or in the beam synthesizer, the two beams undergo a 
superposition. This can be achieved by the rotation matrix like the one 
given in Eq.()37|) [3T]. For instance, if the angle 6 is 90°, the rotation matrix 
takes the form 

If this matrix is applied to the column vector of Eq. , the result is 



(53) 



V2 + 



(54) 



The upper and lower components show the interferences with negative and 
positive signs respectively. 

We have shown previously j2] that the four-by-four transformation ma- 
trices applicable to the Stokes parameters are like Lorentz-transformation 
matrices applicable to the space-time Minkowskian vector {ct,z,x,y). This 
allows us to study space-time symmetries in terms of the Stokes parameters 
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which are apphcable to interferometers. Let us first see how the rotation 
matrix of Eq. (j37|) is translated into the four- by-four formahsm. In this case, 



a = 6 = cos(^/2), 7 = -/? = sin(^/2). 
The corresponding four-by-four matrix takes the form ^ 



R{e) 



(55) 



/I 








°\ 





COS 6* 


— sin 6 








sin 6' 


cos 6 





Vo 








1/ 



(56) 



Let us next see how the phase-shift matrix of Eq. (j5(jp is translated into 
this four-dimensional space. For this two-by-two matrix, 



a 



-i(j}/2 



/5 = 7 = 0, 6 = e^^/2_ (57) 
For these values, the four-by-four transformation matrix takes the form ^ 

\ 



P(0) 



/I 

1 

COS ( 

V sin ( 






— sin I 

cos (h 



(58) 



For the squeeze matrix of Eq. , 



a 



7 



0, 



6 



--0/2 



(59) 



As a consequence, its four-by-four equivalent is 

/ cosh r] 
sinhi] 


V 



sinh 7] \ 

cosh 1] 

1 

Ij 







(60) 



If the above matrices are applied to the four-dimensional Minkowskian space 
of {ct,z,x,y), the above squeeze matrix will perform a Lorentz boost along 
the z or 5*1 axis with 5*0 as the time variable. The rotation matrix of Eq. ()56p 
will perform a rotation around the y or 5*3 axis, while the phase shifter of 
Eq. lj^Hj) performs a rotation around the z or the 5*1 axis. Matrix multipli- 
cations with R{6) and P{(f)) lead to the three-parameter group of rotation 
matrices applicable to the three-dimensional space of {Si, S2, S^,). 
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The phase shifter P(0) of Eq. (j58|) commutes with the squeeze matrix of 
Eq.(jnni), but the rotation matrix R{9) does not. This aspect of matrix algebra 
leads to many interesting mathematical identities which can be tested in 
laboratories. One of the interesting cases is that we can produce a rotation 
by performing three squeezes jl]. Another interesting case is a combination 
of squeeze and rotation matrices which will lead to a triangular matrix with 
unit diagonal elements. This aspect is known as the Iwasawa decomposition 
and is discussed in detail in Ref. jl]. 

7 Density Matrices and Their Little Groups 

According to the definition of the density matrix ^j, the coherency matrix 
of Eq. (j45|) is also the density matrix. In this section, we shall discuss the 
cohrency matrix as the density matrix. 

Under the influence of the G transformation given in Eq. ()46p . this co- 
herency matrix is transformed as 



where C and G are the density matrix and the transformation matrix given 
in Eq. fl45j) and Eq. ()4(ij) respectively. According to the basic property of the 
Lorentz group, these transformations do not change the determinant of the 
density matrix C. Transformations which do not change the determinant are 
called unimodular transformations. 

As we shall see in this section, the determinant for pure states is zero, 
while that of mixed states does not vanish. Is there then a transformation 
matrix which will change this determinant within the Lorentz group. The 
answer is No. This is the basic issue we would like to address in this section. 

If the phase difference between the two waves remains intact, the system 
is said to in a pure state, and the density matrix can be brought to the form 



through the transformation of Eq. with a suitable choice of the G matrix. 





(61) 




(62) 
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For the pure state, the Stokes four-vector takes the form 



1 




(63) 



In order to study the symmetry properties of the density matrix, let us 
ask the following question. Is there a group of transformation matrices which 
will leave the above density matrix invariant? In answering this question, it is 
more convenient to use the Stokes four- vector. The column vector of Eq. (j63|) 
is invariant under the operation of the phase shifter P(0) of Eq.()58p. In 
addition, it is invariant under the following two matrices: 



Fiiu) 



F2(V) 



(l+u'/2 




u 








I u 





u 


—u 


1 





V 








1/ 


(l+v''/2 


-v^/2 





.\ 


v^/2 


1 - v^/2 





V 








1 





V u 


—V 





1/ 



(64) 



These mathematical expressions were first discovered by Wigner in 1939 
in connection with the internal space-time symmetries of relativistic parti- 
cles. They went through a stormy history, but it is gratifying to note that 
they serve a useful purpose for studying interferometers where each matrix 
corresponds to an operation which can be performed in laboratories. 

The F\ and matrices commute with each other, and the multiplication 
of these leads to the form 





[u^ + v^)/2 
iu^ + v^)/2 


u 


u\ 


(M2 + t;2)/2 1- 


u 


V 


U 


—u 


1 





\ V 


—V 





1/ 



(65) 



This matrix contains two parameters. 

Let us go back to the phase-shift matrix of Eq. (j58j) . This matrix also 
leaves the Stokes vector of Eq. (j63j) invariant. If we define the "little group" 
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as the maximal subgroup of the Lorentz group which leaves a Stokes vector 
invariant, the httle group for the Stokes vector of Eq. (jHH|) consists of the 
transformation matrices given in Eq. ()58p and Eq. (jU3j) . 

Next, if the phase relation is completely random, and the first and second 
components have the same amplitude, the density matrix becomes 



1/2 






1/2 



(66) 



Here is the question: Is there a two-by-two matrix which will transform the 
pure-state density matrix of Eq. ()62|) into the impure-state matrix of Eq. ()66|) ? 
The answer within the system of matrices of the form given in Eq. PU)) is No, 
because the determinant of the pure-state density matrix is zero while that 
of the impure-state matrix is 1/4. Is there a way to deal with this problem? 
This problem was addressed in Ref. [31]. In this section, we restrict ourselves 
to the unimodular transformation of Eq. ()61|) which preserves the value of the 
determinant of the density matrix. The Stokes four- vector corresponding to 
the above density matrix is 



\0J 

This vector is invariant under both the rotation matrix of Eq. ()56p and the 
phase shift matrix of Eq.(|3Hj). Repeated applications of these matrices lead 
to a three-parameter group of rotations applicable to the three-dimensional 
space of {Si, S2, S3). 

Not all the impure-state density matrices take the form of Eq. ()66|) . In 
general, if they are brought to a diagonal form, the matrix takes the form 



(67) 



1 + cos X 








1 — COS X 

and the corresponding Stokes four- vector is 

/ cosh fj \ 
sinh?] 





(68) 



v 



with 



1 1 + cos X 

2 1 — cos X 



(69) 



(70) 
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The matrix which transforms Eq. (j67|) to Eq. (j69|) is the squeeze matrix of 
Eq . (|6(J|) . The question then is whether it is possible to transform the pure 
state of Eq. ljUHj) to the impure state of Eq. ljU^ or to that of Eq. (jU7j) . 

In order to see the problem in terms of the two-by-two density matrix, let 
us go back to the pure-state density matrix of Eq.(|U^. Under the rotation 
of Eq.(jS7I), 

/cos(x/2) -sin(x/2)Wl 0Wcos(x/2) sin(x/2) \ . . 

Uin(x/2) cos(x/2) Ao oA-sin(x/2) cos(x/2) ^ ' ^ ^ 

the pure-state density matrix becomes 

1/1 + cosx sinx \ ^^2) 
2 \ sin X 1 — cos xJ ' 

For the present case of two-by-two density matrices, the trace of the 
matrix is one for both pure and impure cases. The trace of the {matrix^ is 
one for the pure state, while it is less than one for impure states. 

The next question is whether there is a two- by-two matrix which will 
eliminate the off-diagonal elements of the above expression that will also 
lead to the expression of Eq. ljUHj) . In order to answer this question, let us 
note that the determinant of the density matrix vanishes for the pure state, 
while it is non-zero for impure states. The Lorentz-like transformations of 
Eq. (j6T| leave the determinant invariant. Thus, it is not possible to transform 
a pure state into an impure state by means of the transformations from the 
six-parameter Lorentz group. Then is it possible to achieve this purpose 
using two-by-two matrices not belonging to this group. We do not know the 
answer to this question. We are thus forced to resort to four-by-four matrices 
applicable to the Stokes four-vector. 



8 Decoherence Effects on the Little Groups 

We are interested in a transformation which will change the density matrix of 
Eq. (jn2I) to that of Eq. (jUUI) . For this purpose, we can use the Stokes four- vector 
consisting of the four elements of the density matrix. The question then is 
whether it is possible to find a transformation matrix which will transform the 
pure-state four- vector of Eq. ljUHj) to the impure-state four- vector of Eq.(jnZj). 

Mathematically, it is more convenient to ask whether the inverse of this 
process is possible: whether it is possible to transform the four-vector of 
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/ cosh 7] 


sinh rj 











/ cosh rj \ 


sinh 1] 


cosh rj 













sinh rj 








1 













V 








l) 






V J 



Eq.(|67|) to that of Eq. (j63|) . This is known in mathematics as the contraction 
of the three-dimensional rotation group into the two-dimensional Euclidean 
group fTTj . Let us apply the squeeze matrix of Eq. (plj) to the four- vector of 
Eq. (jU7j) . This can be written as 



(73) 



After an appropriate normalization, the right-hand side of the above equation 
becomes like the pure-state vector of Eq.(|63p in the limit of large t], as cosh 77 
becomes equal to sinh// in the infinite-77 limit. This transformation is from 
a mixed state to a pure or almost-pure state. Since we are interested in the 
transformation from the pure state of Eq. ljUnj) to the impure state of Eq. (jU7j) . 
we have to consider an inverse of the above equation: 

/ cosh T] 

— sinh 7] 





V 



— sinh 7] 
cosh?7 






1 









\ 


/ cosh 7] \ 








sinh?7 

















/ 


I J 







(74) 



However, the above equation does not start with the pure-state four-vector. 
If we apply the same matrix to the pure state matrix, the result is 



/ cosh T] 

- sinh 7] 


V 



— sinh 7] 
cosh?7 









1 
ij 



1 





1 





(75) 



The resulting four-vector is proportional to the pure-state four-vector and is 
definitely not an impure-state four-vector. 

The inverse of the transformation of Eq. ()73|) is not capable of bringing 
the pure-state vector into an impure-state vector. Let us go back to Eq. (f7S|) . 
it is possible to bring an impure-state into a pure state only in the limit of 
infinite 7]. Otherwise, it is not possible. It is definitely not possible if we take 
into account experimental considerations. 

The story is different for the little groups. Let us start with the rotation 
matrix of Eq. ()56|) . and apply to this matrix the transformation matrix of 
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Eq.(|72D. Then 



/ cosh rj sinh rj \ 












°\ 


sinh?] cosh 7] 







COS 6 


— sin 6 







10 







sin^ 


cos 9 
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ly^ 












1/ 




/ cosh 1] — sinh t] 





^) 










— sinh rj cosh rj 














X 





1 













^ 





l) 









(76) 



If 7] is zero, the above expression becomes the rotation matrix of Eg . (1561) . 
If rj becomes infinite, it becomes the little-group matrix Fi{u) of Eg . (j64j) 
applicable to the pure state of Eg. ()63|) . The details of this calculation for the 
case of Lorentz transformations are given in the 1986 paper by Han et al. ^3] . 
We are then led to the question of whether one little-group transformation 
matrix can be transformed from the other. 

If we carry out the matrix algebra of Eg. (j76p , the result is 



av?w/2 l — v?w/2 
auw —uw 

V 



auw \ 

uw 

'a^)u'^w/2 

IJ 



(77) 



where 



a 



tanhr]. 



u 




w 



1 + (1 -a2)tan2(^/2)' 



(78) 



If a = 0, the above expression becomes the rotation matrix of Eg. (j56|) . If 
a = 1, it becomes the Fi matrix of Eg.()64|). Here we used the parameter a 
instead of r). In terms of this parameter, it is possible to make an analytic 
continuation from the pure state with a = 1 to an impure state with a < 1 
including a = 0. 

On the other hand, we should keep in mind that the determinant of the 
density matrix is zero for the pure state, while it is non-zero for all impure 
states. For a = 1, the determinant vanishes, but it is nonzero and stays the 
same for all non-zero values of a less than one and greater than or egual to 
zero. The analytic expression of Eg. ()78|) hides this singular nature of the 
little group P^ . 
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9 Poincare Sphere as the Representation of 
the Lorentz group 

In Sees. EJ El and |H1 it was noted that the Stokes parameters form a four- 
vector in the Minkowskian space. Thus, it was possible to discuss the density 
matrix in terms of the four-vectors. 

The Poincare sphere was originally constructed from polarization optics. 
Therefore, it is also a representation of the Lorentz group. The Poincare 
sphere for various polarization states have been thoroughly discussed in a 
recent book by Brosseau [U]. 

What is interesting is that the Poincare sphere has two radii. One of them 
is the maximum radius specified by 5*0, and the other radius is the length of 
the space-like components, namely 



is invariant under Lorentz transformations. On the other hand, M = for a 
fully coherent state, while it is non-zero for a partially coherent state. The 
system is totally incoherent if S* = 0. However, the Lorentz group cannot 
handle this decoherence mechanism. 

We observe here that the pure state with M = corresponds to the E{2)- 
like little group for massless particles while it corresponds to the 0(3)-like 
little group for non-zero values of M. The transition of 0(3) to E{2) is 
well known as the group contraction. However, the inverse transformation 
requires further study. 



In analyzing optical rays in para-axial lens optics, we start with the lens 
matrix: 




(79) 



According to the four-vector property of the Stokes parameter. 



(80) 



10 One-lens System 




(81) 



and the translation matrix 




(82) 
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assuming that the beam is propagating along the z direction. 
Then the one-lens system consists of 



\ Z2\( I OWl 
lAo 1 



(83) 



If we perform the matrix multiplication, 

(84) 



1-^2// ^1 + 22-^1^2// 
-1// 1-^1// 

If we assert that the upper-right element be zero, then 

i + i^. (85) 

and the image is focussed, where zi and Z2 are the distance between the 
lens and object and between the lens and image respectively. They are in 
general different, but we shall assume for simplicity that they are the same: 
zi = Z2 = z. We are doing this because this simplicity does not destroy the 
main point of our discussion, and because the case with two different values 
has been dealt with in the literature J2] • Under this assumption, we are left 
with 

1-z/f 2z-z^lf' 
-1// 1-^// 

The diagonal elements of this matrix are dimensionless. In order to make 
the off-diagonal elements dimensionless, we write this matrix as 



(86) 



z \(l-zlf z/f-2\f^ 
\ z/f l~z/f)\Q 



(87) 



Indeed, the matrix in the middle contains dimensionless elements. The neg- 
ative sign in front is purely for convenience. We are then led to study the 
core matrix 

Here, the important point is that the above matrices can be written in 
terms of transformations of the Lorentz group. In the two-by-two matrix 
representation, the Lorentz boost along the z direction takes the form 
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and the rotation around the y axis can be written as 



and the boost along the x axis takes the form 

- sinh(x/2)\ , . 

Uinh(x/2) cosh(x/2)y'- ^ > 

Then the core matrix of Eq. (j88p can be written as 

B{r^)R{cp)B{-T^), (92) 



or 

cos(0/2) -e"'' sin(0/2) \ 
,e^sin((/)/2) cos((/)/2) ) 



(93) 



if 1 < X < 2. If X is greater than 2, the upper-right element of the core is 
positive and it can take the form 

i?(r/)X(x)5(-r/), (94) 



or 

/ cosh(x/2) e-''sinh(x/2)\ , . 

U^" sinh(x/2) cosh(x/2) 7" ^ > 

The expressions of Eq. ()92|) and Eq. ()94|) are a Lorentz boosted rotation 
and a Lorentz-boosted boost matrix along the x direction respectively. These 
expressions play the key role in understanding Wigner's little groups for 
relativistic particles. 

Let us look at their explicit matrix representations given in Eq.()93|) and 
Eq.fl^. The transition from Eq.(|^ to Eq.(|^ requires the upper right 
element going through zero. This can only be achieved through rj going to 
infinity. If we like to keep the lower-left element finite during this process, 
the angle and the boost parameter x have to approach zero. The process 
of approaching the vanishing upper-right element is necessarily a singular 
transformation. This aspect plays the key role in unifying the internal space- 
time symmetries of massive and massless particles. This is like Einstein's 
E = \J {pcY + mgC^ becoming E = pc'm. the limit of large momentum. 

On the other hand, the core matrix of Eq. ()88|) is an analytic function 
of the variable x. Thus, the lens matrix allows a parametrization which 
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allows the transition from massive particle to massless particle analytically. 
The lens optics indeed serves as the analogue computer for this important 
transition in particle physics. 

/,From the mathematical point of view, Eq. ljUHj) and Eq. ()95|) represent 
circular and hyperbolic geometries, respectively. The transition from one 
to the other is not a trivial mathematical procedure. It requires a further 
investigation. 

Let us go back to the core matrix of Eq. ljHHj) . The x parameter does not 
appear to be a parameter of Lorentz transformations. However, the matrix 
can be written in terms of another set of Lorentz transformations. This 
aspect has been discussed in the literature [Si- 



ll Multi-lens Problem 

Let us consider a co-axial system of an arbitrary number of lens. Their focal 
lengths are not necessarily the same, nor are their separations. We are then 
led to consider an arbitrary number of the lens matrix given in Eq.()81|) and 
an arbitrary number of translation matrix of Eq.(|82j). They are multiplied 
like 

Ti Li T2 L2 T3 L3 Tn Ln, (96) 

where N is the number of lenses. 

The easiest way to tackle this problem in to use the Lie-algebra approach. 
Let us start with the generators of the Sp(2) group: 

Since the generators are pure imaginary, the transformation matrices are real. 
On the other hand, the L and T matrices of Eq. (j81|) and Eq.(|82j) are 

generated by 

i\ _ /O 0' 



If we introduce the third matrix 

X,= {1 '_;}, (99) 

all three matrices form a closed set of commutation relations: 

[Xi, X2] = tXs, [Xi, X3] = -^Xi, [X2, X3] = tX2. (100) 
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Thus, these generators also form a closed set of Lie algebra generating real 
two-by-two matrices. What group would this generate? The answer has to 
be 5*^(2). The truth is that the three generators given in Eq. pOOj) can be 
written as linear combinations of the generators of the Sp{2) group given in 
Eq.(jnZ|) 7\. Thus, the Xi matrices given above can also act as the generators 
of the Sp{2) group, and the lens-system matrix given in Eq. (j96p is a three- 
parameter matrix of the form of Eq. (j39|) with real elements. 
The resulting real matrix is written as 



A B 
C D 



(101) 



and is called the ABCD matrix. The question then is how many lenses are 
need to give the most general form of the ABCD matrix [33]. This matrix 
has three independent parameter. 

According to Bargmann ^SGj, this three-parameter matrix can be decom- 
posed into 

>) -sin(a/2)\ Wcos(/?/2) - sin(/5/2) \ 

>) cos(«/2) A e-^/V Uin(/?/2) cos(/?/2) ^ ' ^ ' 

which can be written as the product of one symmetric matrix resulting from 

(103) 



cos(a/2) -sin(a/2) g \ / cos(a/2) sin(a/2) 

sin(a/2) cos(a/2) ) \ Q e"^!'^ )\- sin(a/2) cos(a/2) 

and one rotation matrix: 

cos[(/5 + a)/2] -sin[(/3 + a)/2]\ 
sin[(/3 + a)/2] cos[(/5 + a)/2] j ' 



(104) 



The expression in Eq. (jl03|) is a symmetric matrix, while that of Eq. (jl04j) 
is an orthogonal one. We can then decompose each of these two matrices 
into the lens and translation matrices. The net result is that we do not need 
more than three lenses to describe the lens system consisting of an arbitrary 
number of lenses. The detailed calculations are given in Ref . [7] . 

12 Laser Cavities 

In a laser cavity, the optical ray makes round trips between two mirrors. 
One cycle is therefore equivalent to a two-lens system with two identical 
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lenses and the same distance between the lenses. Let us rewrite the matrix 
corresponding to the one-lens system given in Eq. (j88j) . 



C=f^ ^ ^ ?V (105) 
\ X X — 1 J 

Then one complete cycle consists of C^. 

Here, the cycle starts from the midpoint between the two mirrors [7], 
unlike the traditional approach to this problem where the cycle starts from 
the surface of one of the two mirrors jSZ]- By choosing the midpoint, we can 
eliminate the auxiliary fiat mirror between them needed in the traditional 
approach. This problem was discussed in detail in Ref. [7]. 

For cycles, the expression should be C^^. This calculation, using the 
above expression, will not lead to a manageable form. However, we can resort 
to the expressions of Eq. ()92|) and Eq. ()94|) . Then one cycle consists of 

C = [Biv)R{cjy/2)B{-v)][B{r])R{cl>/2)B{-v)] 

= B{7])R{(l))B{-r]), (106) 

if the upper-right element is negative. If it is positive, the expression should 
be 

C = [Biv)Xix/2)Bi-r^)][Biv)XixmB{-v)] 

= B{r^)X{x)B{-r,). (107) 

If these expressions are repeated times, 

C^^ = B{7])R{N(f))B{-r]), (108) 

if the upper-right element is negative. If it is positive, the expression should 
be 

C'"^ = Bir])XiNx)Bi-7^). (109) 

As A^ becomes large, cosh(A^x) and sinh(A^x) become very large, the beam 
deviates from the laser cavity. Thus, we have to restrict ourselves to the case 
given in Eq. ()108p . 

The core of the expression of Eq. ()108|) is the rotation matrix 

RiN<P) = [R{<P)f. (110) 
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This means that one complete cycle in the cavity corresponds to the rotation 
matrix R{4>). The rotation continues as the beam continues to repeat the 
cycle. 

Let us go back to Eq. ()106|) . The expression corresponds to a Lorentz 
boosted rotation, or bringing a moving particle to its rest frame, rotate, 
and boost back to the original momentum. The rotation associated with 
the momentum-preserving transformation is called the Wigner's little-group 
rotation, which is related to the Wigner rotation commonly mentioned in the 
literature i7||. 



13 Multi-layer Optics 



;iii) 



The most efficient way to study multi-layer optics is to use the S-matrix 
formalism jTUj. This formalism is also based on two-by-two matrices, and we 
can write 

.^J \c dJ\o 

where ipi,ip2 and ip3 are the incoming, reflected and transmitted beams. Here 
we use the matrix notion ABCD for the two-by-two 5* matrix. 

We consider a system of two different optical layers. For convenience, we 
start from the boundary from medium 2 to medium 1. We can write the 
boundary matrix as 



R^r.^ _ /^cosh(r//2) sinh(r7/2)\ , . 

Umh(V2) cosh(r7/2)y' ' ^ ' 

taking into account both the transmission and reflection of the beam. The 
parameter r] is of course determined by the transmission and reflection co- 
efficient at the surface. This problem was studied in depth by Monzon et 
al. jHHl; and their research on this subject is still continuing. 

As the beam goes through the medium 1, the beam undergoes the phase 
shift represented by the matrix 

Wi)=(\ (113) 
When the wave hits the surface of the second medium, the boundary matrix 



IS 



f>(-r,\-( cosh(?7/2) -sinh(r//2)\ , . 

^ ^' l-sinh(r7/2) cosh(r7/2) ^ ' 
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which is the inverse of the matrix given in Eg . (11121) . Within the second 
medium, we write the phase-shift matrix as 

/ -j02/2 Q \ 

m)=('o (115) 

Then, when the wave hits the first medium from the second, we have to go 
back to Eq. (jll2|) . Thus, the one cycle consists of 

cosh(r7/2) sinh(r//2) \ /e-^<^i/2 g 
sinh{r]/2) cosh(r7/2)j V e'^^/^ 



/ cosh(V2) -sinh(V2)\ /e-*^/2 u n . 
^-sinh(r//2) cosh(V2) A p^M^ I ■ l-^^^J 



The above matrices contain complex numbers. However, it is possible to 
transform simultaneously 

/cosh(r7/2) sinh(r//2)\ , . 

Vsinh(r//2) cosh(V2) / ^ ' 

(118) 



to 



/exp(V2) 
V exp (-r//2) 



and transform 



to 



g-i0,/2 Q 

Q g«0i/2 



(119) 
(120) 



'cos(0i/2) -sin(0j/2) 
,sin(0i/2) cos(0i/2) 

using a conjugate transformation, as is shown in Appendix O It is also 
possible to transform these expressions back to their original forms. This 
transformation property has been discussed in detail in Ref . JT] , and also in 
Appendix IHl of the present paper. 

As a consequence, the matrix of Eq. ()116|) becomes 

/cos(0i/2) — sin(0i/2) 
e-''/2 J [ sin(0i/2) cos(0i/2) 



e-''/^ \ /cos(02/2) -sin(02/2)\ 
e''/2 J [ sin(02/2) cos(02/2) ) 



;i2i) 
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In the above expression, the first three matrices are of the same mathematical 
form as that of the core matrix for the one- lens system given in Eq. (j92|) . The 
fourth matrix is an additional rotation matrix. This makes the mathematics 
of repetition more complicated, but this has been done |12j . 
As a consequence the net result becomes 

B{ij)R{Na)B{-ij), (122) 

or 

B{fi)X{NOBi-fi), (123) 

where the parameters fi, a and ^ are to be determined from the input pa- 
rameters ?7, 01 and 02- Detailed calculations are given in Ref. [T^ . 

It is interesting to note that the Lorentz group can serve as a computa- 
tional device also in multi-layer optics. 

14 Concluding Remarks 

We have seen in this report that the Lorentz group provides convenient cal- 
culational tools in many branches of ray optics. The reason is that ray optics 
is largely based on two-by-two matrices. These matrices also constitute the 
group SL{2, c) which serves as the universal covering group of the Lorentz 
group. 

The optical instruments discussed in this report are the fundamental com- 
ponents in optical circuits. In the world of electronics, electric circuits form 
the fabric of the system. In the future high-technology world, optical compo- 
nents will hold the key to technological advances. Indeed, the Lorentz group 
is the fundamental language for the new world. 

It is by now well known that the Lorentz group is the basic language 
for quantum optics. Coherent and squeezed states are representations of 
the Lorentz group. It is challenging to see how the Lorentz nature of the 
above-mentioned optical components will manifest itself in quantum world. 

The Lorentz group was introduced to physics by Einstein and Wigner 
to understand the space-time symmetries of relativistic particles and the 
covariant world of electromagnetic fields. It is gratifying to note that the 
Lorentz group can serve as the language common both to particle physics 
and optical sciences. 
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Appendix 



A Lorentz-squeezed Harmonic Oscillators 

The purpose of this Appendix is to expand the third row of Table [T] using 
the covariant formahsm of harmonic oscillators. We explain first what this 
formalism does for covariant picture of relativistic extended particles. We 
then point out that this oscillator formalism can serve as the basic language 
for two-mode coherent states |2Hl5 namely squeezed states of light. 

The concept of localized probability distribution is the backbone of the 
present formulation of quantum mechanics. Of course, we would like to have 
more deterministic form of dynamics, and efforts have been and are still 
being made along this direction. One of the most serious problems with this 
probabilistic interpretation is whether this concept of probability is consistent 
with special relativity. 

In a given Lorentz frame, we know how to do quantum mechanics with 
a localized probability distribution. How would this distribution look to an 
observer in a different Lorentz frame? 

• Would the probability distribution appear the same to this observer? 

• If different, how is the probability distribution distorted? 

• Is the total probability conserved? 

• The Lorentz boost mixes up the spatial coordinate with the time vari- 
able. What role does the time-separation variable play in defining the 
boundary condition for localization and the probability distribution? 

We can answer some or all of the above questions only if we construct co- 
variant wave functions, namely wave functions which can be Lorentz-boosted. 
It is easy to construct these wave functions if we know the answer. In the 
initial development of quantum mechanics, the harmonic oscillator played 
the pivotal role. Thus, it is clear to us that if there is a wave function which 
can be Lorentz-boosted, this has to be the harmonic oscillator wave func- 
tion. Until we construct wave functions which can be Lorentz-transformed, 
we cannot say that quantum mechanics is consistent with relativity. Indeed, 
we should examine this problem before attempting to invent more definitive 
quantum mechanics. 
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Since the hadron, in the quark model, is a bound-state of quarks PU] , 
Feynman, Kishnger, and Ravndal, in 1971 [lU]; raised the following question 
in connection with the quark model for hadrons: the hadronic spectrum 
can indeed be explained in terms of the degeneracy of three-dimensional 
harmonic oscillator wave functions in the hadronic rest frame; however, what 
happens when the hadron moves? Indeed, Feynman et al. wrote a Lorentz- 
invariant differential equation whose solutions can become non-relativistic 
wave functions if the time-separation variable can be ignored. 

This Lorentz-invariant differential equation is a four-dimensional partial 
differential equation, with many different solutions depending on separation 
of variables and boundary conditions. There is a set of normalizable solutions 
which can serve as a representation space for Wigner's little group for massive 
particles [U [T7j. We can start with this set of solutions and give physical 
interpretations, especially to the time-separation variable. 

Finally, is the covariance of the oscillator wave function consistent with 
what we observe in the real world. Here again, Feynman plays the key 
role. While the quark model can be fit into the oscillator scheme in the 
hadronic rest frame, Feynman in 1969 came up with the idea of partons |3T] . 
According to Feynman's parton model, the hadron consists of an infinite 
number of partons when it moves with a velocity close to that of light. Quarks 
and partons are believed to be the same particles, but their properties are 
totally different. While the quarks inside the hadron interact coherently with 
external signals, partons interact incoherently. If the partons are Lorentz- 
boosted quarks, does the Lorentz boost destroy the coherence? 

Before 1964 [Sn], the hydrogen atom was used for illustrating bound 
states. These days, we use hadrons which are bound states of quarks. Let 
us use the simplest hadron consisting of two quarks bound together with an 
attractive force, and consider their space-time positions Xa and Xf,, and use 
the variables 



The four-vector X specifies where the hadron is located in space and time, 
while the variable x measures the space-time separation between the quarks. 
According to Einstein, this space-time separation contains a time-like com- 
ponent which actively participates as can be seen from 





(124) 




(125) 
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when the hadron is boosted along the z direction. In terms of the hght-cone 
variables defined as ^| 



u = 



V = {z-t)/V2, 



(126) 



the boost transformation of Eq. ()125|) takes the form 



u' = e'^u, 



V. 



(127) 



The u variable becomes expanded while the v variable becomes contracted. 

Does this time-separation variable exist when the hadron is at rest? Yes, 
according to Einstein. In the present form of quantum mechanics, we pretend 
not to know anything about this variable. Indeed, this variable belongs to 
Feynman's rest of the universe. In this report, we shall see the role of this 
time-separation variable in the decoherence mechanism. 

Also in the present form of quantum mechanics, there is an uncertainty 
relation between the time and energy variables. However, there are no known 
time-like excitations. Unlike Heisenberg's uncertainty relation applicable to 
position and momentum, the time and energy separation variables are c- 
numbers, and we are not allowed to write down the commutation relation 
between them. 

How does this space-time asymmetry fit into the world of covariance j22j . 
The answer is that Wigner's 0(3)-like little group is not a Lorentz-invariant 
symmetry, but is a covariant symmetry fl]. It has been shown that the time- 
energy uncertainty relation applicable to the time-separation variable fits 
perfectly into the 0(3)-like symmetry of massive relativistic particles |17j . 

The c-number time-energy uncertainty relation allows us to write down a 
time distribution function without excitations [17 . If we use Gaussian forms 
for both space and time distributions, we can start with the expression 



for the ground-state wave function. What do Feynman et al. say about this 
oscillator wave function? 

In their classic 1971 paper |3U], Feynman et al. start with the following 
Lorentz-invariant differential equation. 






(129) 
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This partial differential equation has many different solutions depending on 
the choice of separable variables and boundary conditions. Feynman et al. 
insist on Lorentz-invariant solutions which are not normalizable. On the 
other hand, if we insist on normalization, the ground-state wave function 
takes the form of Eq. ()128|) . It is then possible to construct a representation of 
the Poincare group from the solutions of the above differential equation ^7] . 
If the system is boosted, the wave function becomes 

ij^{z,t) = Q'^'exp|-i(e-2V + e%2)}. (130) 

This wave function becomes Eq . (11281) if r] becomes zero. The transition from 
Eq. (jl28p to Eq. (jl3(J|) is a squeeze transformation. The wave function of 
Eq. ()128|) is distributed within a circular region in the uv plane, and thus in 
the zt plane. On the other hand, the wave function of Eq. ()130|) is distributed 
in an elliptic region with the light-cone axes as the major and minor axes 
respectively. If t] becomes very large, the wave function becomes concen- 
trated along one of the light-cone axes. Indeed, the form given in Eq. (|13U|l 
is a Lorentz- squeezed wave function. This Lorentz-squeeze mechanism is 
illustrated in Fig. [TJ 

There are many different solutions of the Lorentz invariant differential 
equation of Eq. (jl29|) . The solution given in Eq. (jl3Up is not Lorentz invariant 
but is covariant. It is normalizable in the t variable, as well as in the space- 
separation variable z. How can we extract probability interpretation from 
this covariant wave function? This issue has been discussed thoroughly in 
the literature [T7] . 

Another pressing problem in physics is that hadrons, like the proton, 
can be regarded as quantum bound states of quarks when they move slowly. 
However, they appear like collections of partons when they move with velocity 
very close to that of light jH]. Since the quarks and partons have quite 
different properties, it is a challenging problem in modern physics to show 
that they are the same but appear different depending on the observer's 
Lorentz frame. 

It has been shown that this Lorentz-squeezed wave function can explain 
the quark model and the parton model as two different manifestations of one 
covariant entity Thus, the covariant harmonic oscillator can occupy the 
third row of Table ^ 
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At 




Figure 1: Effect of the Lorentz boost on the space-time wave function. The 
circular space-time distribution at the rest frame becomes Lorentz-squeezed 
to become an elhptic distribution. The first version of this figure consists of 
two elhpses in the 1973 paper by Kim and Noz j^Tj. This figure consisting 
of a circle and a ellipse is from the 1978 paper by Kim, Noz, and Oh |43j . 
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The Lorentz-squeezed wave function can be written as 



IT 



i)ri{z, t) = ^ exp <^ -- e^iz - tf + e-'^iz + tf 



(131) 



As was discussed in the hterature for several different purposes [21 El IS]) 
this wave function can be expanded as 



= (tanh|)%fc(2;)0,(t). (132) 



This is an expansion in terms of the oscillations along two different directions. 
If those directions are xi and X2, we can write Eq. p32|) as 



1 ^ / ,7]'^^ 



cosh T] ^ 



E tanh- 0fc(xi)0fc(x2). 



(133) 



The mathematics of harmonic oscillators can be translated into second- 
quantized Fock space. In the Fock space, this expression describes the two- 
photon coherent state discussed first by Yuen in 1976 |2H|- 



B Euler versus Lie Representations 

The group SL{2, c) consists of two-by-two unimodular matrices whose ele- 
ments are complex. There are therefore six independent parameters, and 
thus six generators of the Lie algebra. This group is locally isomorphic to 
the six-parameter Lorentz group or 0(3, 1) applicable to the Minkowskian 
space of three space-like directions and one time-like direction. 

We can construct representations of the group starting from the Lie al- 
gebra given in Sees. El and IHl We can construct the representation by using 
a method similar to what Goldstein did for the three-dimensional rotation 
group in terms of the Euler angles JH]- There are three-generators for the 
rotation group, but Goldstein starts with rotations around the z and x direc- 
tions. Rotations around the y axis and the most general form for the rotation 
matrix can be constructed from repeated applications of those two starting 
matrices. Let us call this type of approach the "Euler construction." 

In constructing the Lorentz group, we observe first that we need rotations 
around two different directions. As for the boost, we need boosts along 
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one given direction since boots along other directions can be achieved by 
rotations. 

There are three basic advantages of this approach. First, the number of 
"starter" matrices is less than the number of generators. For example, we 
need only two starters for the three-parameter rotation group. In our case, 
we started with two matrices for the three-parameter group 5*^(2) and also 
for SU{1,1). Second, each starter matrix takes a simple form and has its 
own physical interpretation. 

The third advantage can be stated in the following way. Repeated ap- 
plications of the starter matrices will lead to a very complicated expression. 
However, the complicated expression can be decomposed into the minimum 
number of starter matrices. For example, this number is three for the three- 
dimensional rotation group. This number is also three for SU{2) and Sp{2). 
We call this the Euler decomposition. The present paper is based on both 
the Euler construction and the Euler decomposition. 

Among the several useful Euler decompositions, the Iwasawa decompo- 
sition plays an important role in the Lorentz group. We have seen in this 
paper what the decomposition does to the two- by-two matrices of Sp{2), but 
it has been an interesting subject since Iwasawa's first publication on this 
subject [21]. It is beyond the scope of this paper to present a historical review 
of the subject. However, we would like to point out that there are areas of 
physics where this important mathematical theorem was totally overlooked. 

For instance, in particle theory, Wigner's little groups dictate the inter- 
nal space-time symmetries of massive and massless particles which are locally 
isomorphic to 0(3) and -£'(2) respectively pp. The little group is the maxi- 
mal subgroup of the Lorentz group whose transformations do not change the 
four- momentum of a given particle [211. The i?(2)-like subgroup for mass- 
less particles is locally isomorphic to the subgroup of SL{2, c) which can 
be started from one of the matrices in Eq. ()149|) and the diagonal matrix of 
Eq. ()113|) . Thus there was an underlying Iwasawa decomposition while the 
£'(2)-like subgroup was decomposed into rotation and boost matrices 
but the authors did not know this. 

In optics, there are two-by-two matrices with one vanishing off-diagonal 
element. It was generally known that this has something to do with the 
Iwasawa effect, but Simon and Mukunda |17] and Han et al. j3] started 
treating the Iwasawa decomposition as the main issue in their papers on 
polarized light. 

In para-axial lens optics, the matrices of the form given in Eq. (jl49|) are 
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the starters, and repeated applications of those two starters will lead to the 
most general form of Sp{2) matrices. It had been a challenging problem 
since 1985 [33] to write the most general two-by-two matrix in lens optics 
in terms the minimum number of those starter matrices. This problem has 
been solved recently [7j , and the central issue in the problem was the Iwasawa 
decomposition. 



C Conjugate Transformations 

The core matrix of Eq. ()88|1 contains the chain of the matrices 

W= (^''^ ° ^ /coshr/ sinhr^X f e"'^ \ 
V e"^ ) \ sinh rj cosh r/ / V e* V ' 

The Lorentz group allows us to simplify this expression under certain condi- 
tions. 

For this purpose, we transform the above expression into a more conve- 
nient form, by taking the conjugate of each of the matrices with 



Then CiWC^^ leads to 

(136) 



cos (f) — sin \ / cosh r] sinh r]\ / cos ^ — sin ^ 
sin cos / V sinh i] cosh i] J \ sin ^ cos ^ 

In this way, we have converted W of Eq. ()134j) into a real matrix, but it is 
not simple enough. 

Let us take another conjugate with 



Then the conjugate C2C1WC1 C2 becomes 

cos — sin \ / e'' ^ \ ( '^'^^ ^ ~ sin ^ 
sin0 COS0 / V e^''/ ysin^ cos^ 



(138) 



The combined effect of C2C1 is 

C = C2C1 = y_p-in/4 p-in/4. ) ? (^39) 
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with 

1 ^g-i7r/4 _gi7r/4 



After muhiphcation, the matrix of Eq. ()138|1 will take the form 

^-{c d)' 

where A, B, C, and D are real numbers. If B and C vanish, this matrix 
will become diagonal, and the problem will become too simple. If, on the 
other hand, only one of these two elements become zero, we will achieve a 
substantial mathematical simplification and will be encouraged to look for 
physical circumstances which will lead to this simplification. 

Let us summarize. We started in this section with the matrix representa- 
tion W given in Eq. ()lH4j) . This form can be transformed into the V matrix 
of Eq. ()lH8j) through the conjugate transformation 

V = CWC-\ (142) 

where C is given in Eq. pifflj) . Conversely, we can recover the W representa- 
tion by 

W = C-^VC. (143) 

For calculational purposes, the V representation is much easier because we 
are dealing with real numbers. On the other hand, the W representation is 
of the form for the S-matrix we intend to compute. It is gratifying to see 
that they are equivalent. 

Let us go back to Eq. ()lH8|l and consider the case where the angles and 
^ satisfy the following constraints. 

+ e = 2^, 0-e = 7r/2, (144) 

thus 

= ^ + 7r/4, ^ = ^-7r/4. (145) 

Then in terms of ^, we can reduce the matrix of Eq. p38|l to the form 

(cosh rf) cos(2^) sinh — (cosh r]) sin(2^) 

sinh ?7 + (cosh r]) sin(2^) (cosh r]) cos(2^) 



(146) 



Thus the matrix takes a surprisingly simple form if the parameters Q and r\ 
satisfy the constraint 

sinh?7 = (coshr^) sin(2^). (147) 
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Then the matrix becomes 



1 
2 sinh rj 1 



148) 



This aspect of the Lorentz group is known as the Iwasawa decomposition j2H] , 
and has been discussed in the optics hterature [H [TTJ H?] . 

Matrices of this form are not so strange in optics. In para-axial lens 
optics, the translation and lens matrices are written as 

; ;). -'^ a ?). (m 

respectively. These matrices have the following interesting mathematical 
property [Sj, 

1 Ui\fl U2\ fl Ui+U2\ ^^^Q^ 



1 y vo 1 y vo i 

and 

' ?)f' ?) = f ' ?V (151) 

We note that the multiplication is commutative, and the parameter becomes 
additive. These matrices convert multiplication into addition, as logarithmic 
functions do. 
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